Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



o 
o 

(N 
o 

Q 

oo 



o 
o 

u 

Oh' 



Kinetic Energy, Condensation Energy, Optical Sum Rule 
and Pairing Mechanism in High-Tc Cuprates 
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The mechanism of high-Tc superconductivity is investigated with interests on the micro- 
scopic aspects of the condensation energy. The theoretical analysis is performed on the basis 
of the FLEX approximation which is a microscopic description of the spin-fluctuation-induced- 
superconductivity. Most of phase transitions in strongly correlated electron system arise from 
the correlation energy which is copmetitive to the kinetic energy. However, we show that the 
kinetic energy cooperatively induces the superconductivity in the underdoped region. This un- 
usual decrease of kinetic energy below Tc is induced by the feedback effect. The feedback effect 
induces the magnetic resonance mode as well as the kink in the electronic dispersion, and alters 
the properties of quasi-particles, such as mass renormalization and lifetime. The crossover from 
BCS behavior to this unusual behavior occurs for hole dopings. On the other hand, the decrease 
of kinetic energy below Tc does not occur in the electron-doped region. We discuss the relation 
to the recent obserbation of the violation of optical sum rule. 

KEYWORDS: High-Tc superconductivity; spin fluctuation; feedback effect; kinetic energy; condensation 
energy; optical sum rule 
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1. Introduction 

The mechanism of high-Tc superconductivity in 
cuprate materials^ has been one of the most appeahng 
subject in the condensed matter physics over the last cou- 
ple of decades. Through intensive studies from theoreti- 
cal and experimental researchers, the "magnetic mecha- 
nism" is believed most predominantly. 

The "magnetic mechanism" is represented by the spin 
fluctuation theory which takes into account the inter- 
action between quasi-particles exchanging the spin fluc- 
tuations.^^ ^ Among the microscopic descriptions beyond 
the original phenomenology, the fluctuation-exchange 
(FLEX) approximation is adopted most widely. ^"^^ The 
qualitative validity of this theory is highly expected from 
weak to intermediate coupling region. -'^'^ For example, 
the spin fluctuation theory is robust for the vertex cor- 
rections arising from the non-RPA tcrms^^' as well 
as those from the multiple spin fluctuation exchange 
terms. On the other hand, in the expansion from 
the strong coupling limit, the resonating valence bond 
(RVB) theory^^"^^ and some numerical methods^^"^^ 
have concluded the d^^ _y2-wave superconductivity where 
the super- exchange interaction plays an essential role. 

In the spin fluctuation theory, Eliashberg equation 
is used for an analysis of the superconductivity. The 
Eliashberg equation provides a clear understanding for 
the mechanism of superconductivity on the basis of the 
BCS picture. Then, the attractive interaction leading to 
the Cooper pairing is represented by the irreducible four 
point vertex. The momentum dependence of this vertex 
induces the non-s-wave superconductivity. In the spin 
fluctuation theory for high-Tc cuprates, the irreducible 
four point vertex is described by the anti-ferromagnetic 
spin fluctuation and the attractive interaction is most 
effective in the d-wave channel. 

The goal of this paper is to investigate the mechanism 



of high-Tc superconductivity from another point of view. 
We study how the energy is gained below Tc owing to the 
superconductivity. Of course, the ground state is deter- 
mined as a result of the energetic optimization. The un- 
derstanding from the energetics will be complementary 
to the analysis of the interaction leading to the pairing. 
Some interesting aspects are clarified from this point of 
view. 

This study is partly motivated by the theoretical pro- 
posal for "kinetic energy driven pairing" and by 
the recent experimental supports for this proposal.''^' 
In the conventional BCS theory, the kinetic energy in- 
creases owing to the superconductivity. This increase 
is slightly over-compensated by the decrease of corre- 
lation energy. Contrary to the BCS theory, the "ki- 
netic energy driven pairing" attributes the mechanism 
of superconductivity to the gain of kinetic energy. This 
mechanism has been considered to be highly uncon- 
ventional and the discrepancy to the spin-fluctuation- 
induced-superconductivity has been noted. ^^"'^'^ On the 
other hand, the consistency to the RVB theory has been 
discussed, where the superconducting transition is trig- 
gered by the coherence of charge carriers.^^''^^''^^ As 
for numerical studies, the dynamical cluster approxima- 
tion'^^ and variational Monte Carlo simulation"^^ for Hub- 
bard model have shown a decrease of kinetic energy ow- 
ing to the superconductivity, and then implications for 
the RVB state have been noted. 

In this paper, we study these problems on the ba- 
sis of the microscopic and strong coupling theory on 
the spin-fluctuation-induced-superconductivity, namely 
the FLEX approximation. This subject has been in- 
vestigated by several phenomenological theories assum- 
ing the non- Fermi liquid normal state, ^^'^^ spin-Fermion 
coupling, superconducting phase fluctuation,"*^ and 
electron-phonon coupling."*^ In contrast to these theo- 
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ries, the FLEX approximation is a "conserving approxi- 
mation" which is highly suitable for a discussion of ther- 
modynamic properties. In the conserving approximation 
formulated by Luttinger and Ward,^'^'^'' all quantities 
are self-consistently derived from the differential of ther- 
modynamic potential without any phenomenological as- 
sumption. Unphysical results inherent in the phenomeno- 
logical theory are considerably excluded in the micro- 
scopic theory adopted here. It should be stressed that 
highly careful treatment is needed for thermodynamic 
properties rather than for magnetic or single-particle 
properties. This is partly because the condensation en- 
ergy of superconductivity is much smaller than the en- 
ergy scale of electrons. For example, the condensation 
energy is in the order of O.lmeV, while the band width 
is in the order of leV. 

Note that kinetic energy along c-axis has attracted 
much interests in the early stage because it is related 
to the 'interlayer tunneling mechanism' (ILT) proposed 
by Anderson.^^ Optical measurements have supported 
the decrease of c-axis kinetic energy. However, it has 
been shown that the gain of c-axis kinetic energy is much 
smaller than the condensation energy,^^"^^ and therefore 
this subject is not essential for the mechanism of super- 
conductivity. In this paper we focus our attention on the 
kinetic energy along the plane. 

In §2, we formulate a conserving approximation in 
the superconducting state and provide the expressions of 
FLEX approximation. Results on the kinetic energy are 
shown in §3.1. We show that the kinetic energy decreases 
below Tc in the under-doped region while it increases like 
BCS theory in the over-doped region and in the electron- 
doped region. It will be stressed that the concepts 
of "spin-fluctuation-induced-superconductivity" and "ki- 
netic energy driven pairing" are not incompatible. The 
relation between the kinetic energy and the optical sum 
rule is discussed in §3.2. In §3.3, wc discuss thermody- 
namic properties in more details. Then, wc propose an- 
other interpretation of the condensation energy by con- 
sidering the free energy arising from the spin fluctuation. 
Some discussions are given in §4. 

2. Thermodynamic Property and FLEX Ap- 
proximation below Tc 

In this paper, we analyze the two-dimensional Hub- 
bard model which is expressed as, 

H = J2^{k)c^l^^C)^^ + UY,n'i^nii. (1) 

We consider the square lattice and choose the following 
tight-binding dispersion, 

e{k) = — 2f(cos kx + cos ky) + At' cos fcx cos fey. (2) 

In the following, the unit of energy is chosen as 2t = 
1. The next nearest neighbor hopping t' is necessary 
and sufficient to reproduce the Fermi surface of high- 
Tc cuprates. The typical value is estimated to be t'/t = 
0.1 ~ 0.4. Qualitative results in this paper are not al- 
tered by this value. We fix t'/t = 0.25 in the hole-doped 
region and t'/t = 0.35 in the electron-doped region, re- 
spectively. The concentration of hole doping is expressed 



as (5 = 1 — n where n is the density of electrons per sites. 

In the superconducting state, statistical quantum field 
theory is described by normal and anomalous Green 
functions, G{k) and F{k). The Dyson-Gorkov equation 
describes the Green functions through the normal and 
anomalous self-energies, which are denoted as S'^(A;) and 
A(fc), respectively. 



G{k) 
Ft (A;) 



m 

-G{-k) 



A*(fc) 



A(fc) I . . 



Here, G^°-* (fc) is the Green function in the non-interacting 
case, 



(4) 



where /i is the chemical potential. The superconduct- 
ing gap A(A;) is obtained by the anomalous self- 
energy as A(fe) = 2;(A;)|A(fc)| where z{k)~^ = 1 — 

dRcE^{k,uj)/duj\^=a- 

In order to discuss the thermodynamic quantities, we 
first formulate a general expression for the thermody- 
namic potential in the superconducting state and derive 
the self-energy, momentum distribution function and ki- 
netic energy on the basis of the functional derivatives. In 
the following, we describe the formulation in case of the 
spin singlet pairing. 

The conserving form of the thermodynamic potential 
in the normal state was formulated by Luttinger and 
Ward^^ and developed by Baym and KadanofF.^^ Then, 
the self-energy is obtained by the functional derivative of 
generating fimction $ as = d^[Ga]/SGcr{k). The 

thermodynamic potential is obtained by the generating 
function and self-energy as. 



n{T,ii) = no{T,ii) 



yVfiogi^^M^^^^} 



+G^{k)^l{k)] + ^G„]. 



(5) 



Here, f2o(T,M) = -2TJ2^log[l + exp{-/3(e(fc) - n)}] 
is the thermodynamic potential in the non-interacting 
case. Although we have formally introduced the index 
of spin a, indeed, G^(fc) = G{k) and S^(A;) = 
since we consider the paramagnetic state or spin singlet 
superconducting state. 

It is straightforward to generalize this formulation 
to the superconducting state. We obtain the normal 
and anomalous self-energies from the generating function 
^[G„,F,F^ as, 

Kik) = 77^, (6) 



A(fc) 



A*(fc) = - 



SG,{k)' 
~SF^k)' 



SF{k) ■ 



(7) 



(8) 



Note that cq. (7) (cquivalcntly eq. (8)) is a self-consistent 
equation determining the second order phase transition. 
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The linearized version of eq. (7) has been used in order to 
determine the superconducting instabihty.^'' We obtain 
the general expression of thermodynamic potential as, 



(9) 



0F = -^[l0g{ 

k 



} 



I -G(0)(fc)-i|2 
+ Y,GAk)K{k) - F{k)A*{k) - F^{k)A{k)UlO) 



(11) 



The derivation of cqs. (9-11) is summarized in Appendix. 
According to eqs. (6-11), the variational conditions with 
respect to the self-energy are satisfied as, 

5n _ set 

S^JJcj ~ SA{k) ~ ^ ' 

Therefore, the self-energies obtained by eqs. (6-8) pro- 
vide a stationary value of thermodynamic potential. Ac- 
cording to the thermodynamics, we obtain the number 
density as 

n=-— = 2^n(fc), (13) 
k 

where the momentum distribution function is also ob- 
tained by the functional derivatives as. 



By performing the functional derivatives, eq. (14) is re- 
duced to the usual definition of n{k), 

n{k) = G{k)e''^-^ 

= ^[G(fc)-G(°)(fc)]+/(£(fc)-M), (15) 

where we have eliminated the ultra-violet divergence by 
subtracting the Fermi distribution function f{e(k) — jj). 
Finally, the kinetic energy is obtained as. 



E^ = Y,e{k) 
k 



sn 

Se{k) 



2^£(fe)n(fc). 
k 



(16) 



When these relations are self-consistently satisfied in 

an approximation, the approximation is classified into 
the "conserving approximation" .^"^ The FLEX approxi- 
mation is one of them. In the following, we fix the number 
density n instead of the chemical potential /i. Therefore, 
the free energy F{T, n) = f2 -t- /un is a more convenient 
quantity. 

Hereafter, the superscripts ^ and ^ denote the siiper- 
conducting and normal states, respectively. The former 
is defined by the stationary solution with finite value 
of A(fc) and the latter is defined by the solution with 
A(fc) = 0. Both solutions satisfy the variational condi- 
tions eq. (12), but only the former satisfies the stable 
condition below T^. The difference between normal and 
superconducting states are denoted as 6 A = — A?. 
For example, the condensation energy is expressed as 
5F = F^ - F^. 



Before closing the general formulation, we note some 
analytical expressions for the parameter dependence of 
the free energy. First, the number dependence of the free 
energy is given by the chemical potential dF/dn = /i. 
Therefore, the number dependence of the condensation 
energy is obtained as, 

dSF „ « 



dn 



(17) 



Second, the t/-dependence of the free energy is given by 
the running coupling constant formula which is expressed 
as follows. 



dU 2U 

-F{k)A*{k) - Ft(fc)A(fc)]. (18) 

The [/-dependence of the condensation energy is simply 
obtained by the subtraction. These expressions are con- 
venient to understand the qualitative behaviors of the 
condensation energy (see §3.3). 

The formulation of FLEX approximation has been 
given in literatures.^' The extension to the supercon- 
ducting state is straightforward. Indeed, sonic; authors 
have investigated the properties in the superconducting 
state by using the FLEX approximation. For instance, 
the temperature dependence of superconducting gap,^'* 
magnetic and single-particle properties^' ^'^^ have been 
discussed. In this paper, we analyze the thermodynamic 
quantities and their relation to the optical sum rule. 

The generating function ^[Ga,F, F'^] is obtained in the 
FLEX approximation as. 



\uHxl{ci)' + Xl(ci)') + U{lx'M - \xlm- (19) 



4^ ...... . xl{qf) + U{-xl{^, 2' 

Here, we have denoted the irreducible spin and charge 
susceptibilities as. 



X, 



k 



[G{k + q)G{k)±F{k + q)F{k)]. (20) 



Note that we have ignored the first order terms in the 
generating function since their roles are trivial and do 
not affect the following discussions. 

We obtain the self-energy from eqs. (6-7) as, 

H^ik) = Y,VM)G{k - q), (21) 

<j 

A{k) = -Y,yM)F{k-q), (22) 



where 14 (g) and Va.{q) are expressed as, 

K(g) = U^ilxM + \xoii) - \{x'M+xlm], (23) 

K(5) = U^[\xM - \xM - \{X%) - X°(9)}]. (24) 
We have introduced the spin and charge susceptibilities 
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obtained by the generalized RPA as, 

As.C V 



Xs,c(«) 



(25) 



The normal vertex Vi-i{q) and anomalous vertex Va((?) are 
represented by the irreducible four point vertex in the 
particle-hole channel and in the particle-particle channel, 
respectively. 

The superconducting transition from the normal state 
is determined by the appearance of non-trivial solution 
of eq. (22). In analogy with the gap equation in the 
BCS theory, we often denote V^{q) as effective interac- 
tion leading to the pairing. Generally speaking, the mo- 
mentum dependence of effective interaction results in the 
attractive interaction in a non-s-wave channel.^'* In the 
FLEX approximation, the effective interaction V^{q) is 
dominated by the spin fluctuation whose momentum de- 
pendence is favorable for the dx2_y2-wave superconduc- 
tivity. This is the ordinary understanding on the spin- 
fluctuation- induced superconductivity. 

In this paper, we obtain another insights on the spin- 
fluctuation-induced superconductivity which are given 
by the analysis of energetics. The free energy is described 
as F = E]i + E„ — TS where i?cr =< Uni^^rn^i > is the 
correlation energy and S is the entropy. At T = 0, the 
condensation energy is obtained by the kinetic energy 
and correlation energy as 5F = 5E\^ + 6Ecr- In §3.1, we 
show that the kinetic energy increases the condensation 
energy cooperatively with the correlation energy. This 
is in sharp contrast with the weak coupling BCS the- 
ory where the kinetic energy remarkably decreases the 
condensation energy. 

Note that the FLEX approximation provides a rea- 
sonable value of Tc not only in the hole-doped region but 
also in the electron-doped region. Tc and doping region 
with superconducting order is very small in the electron- 
doped region^^"^^ owing to the small DOS and the lo- 
calized character of spin fluctuation in the momentum 
space. 

It should be noticed that the FLEX approximation 
does not explain the pseudogap phenomena in the nor- 
mal state of under-doped region. The superconducting 
fluctuation should be included to explain the pseudogap 
phenomena in this framework.^'^'^^'^^ However, it has 
been shown that the effects of superconducting fluctu- 
ation on the electronic state is rapidly suppressed below 
Tc.^* Therefore, we believe that the FLEX approxima- 
tion is appropriate for a description of superconducting 
state. We note that some interesting phenomena, such 
as magnetic resonance peak^^ and kink in the electronic 
dispersion, are well explained within the FLEX approx- 
imation. 

We have used the notations J2k ~ ^/^Sc^ fe ^^'^ 
Eg = T/NJ2Q^,q where w„ = (2n+ l)7rr, f]„ =2mTT, 
T is the temperature and N is the number of sites. 
The unit ft = c = /jb = 1 is used through this paper. 
Eqs. (20), (21) and (22) are estimated by using the fast 
Fourier transformation (FFT). 



3. Results 

Before showing the results, we note some cares in- 
volved in the numerical calculation because a careful 
treatment is highly needed for an estimation of thermo- 
dynamic quantities. This is partly because the electronic 
states far below Fermi level essentially contribute to the 
free energy, and partly because the condensation energy 
is a very small value in the order of 10~* ~ O.lmeV. 
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Fig. 1. (a) Temperature dependence of free energy in the normal 
state calculated by fixing the cut-off frequency. The solid line is 
a fitting curve with use of the function = a + bT'^ + cT'^ log T. 
The inset shows the same quantity with use of fixed Nf = 
1024, 2048 and 4096. (b) Temperature dependence of SF for var- 
ious N{. We choose the parameters as S = 0.1 and U/t = 4.2. 



We divide the first Brillouin zone into N x N and 
take Nf Matsubara frequency. The numerical inaccu- 
racy mainly arises from the cut-off of Matsubara fre- 
quency. It should be noticed that we have introduced 
expressions without ultra-violet divergence in eqs. (9) 
and (15). This procedure remarkably improves the nu- 
merical accuracy. However, further care is needed for 
a temperature dependence of thermodynamic quanti- 
ties. If we fix the number of Matsubara frequency Nf, 
the cut-off of frequency depends on the temperature as 
LOc ~ {Nf — 1)ttT. This induces an artificial temper- 
ature dependence which may smear the intrinsic tem- 
perature dependence. This diflaculty is very serious for 
the free energy, as is shown in Fig. 1(a). The main fig- 
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ure shows the temperature dependence of free energy 
calculated with the cut-off frequency fixed to be 
u)c = 38.6 ~ WW. Then, the free energy is well fitted 
by the function = a + hT"^ + cT^ logT, which is con- 
sistent with the nearly anti-ferromagnetic Fermi liquid 
state.^'^ On the other hand, the free energy obtained 
by fixing iVf shows much larger temperature dependence 
(see inset in Fig. 1(a)) indicating the violation of ther- 
modynamic third law. Thus, we have to fix the cut-off 
frequency instead of Nf in order to obtain appropriate 
results. Unfortunately, owing to the computational con- 
straint arising from the FFT, it is troublesome to fix the 
cut-off frequency u)c = {Nf — 1)'jtT for various tempera- 
tures. 

However, this difficulty does not matter for the differ- 
ences between normal and superconducting states. This 
is because the superconductivity affects the low energy 
states while high energy states are not sensitive to the su- 
perconductivity. For instance, wc show the temperature 
dependence of SF in Fig. 1(b). It is clearly shown that 
SF depends on A^f only slightly, and the fixed- A^f calcula- 
tion is valid. Note that TVf-dependence of the calculated 
free energy is still in the order of 10~^ at T = 0.005. 
However, the difference 5F is estimated to an accuracy 
of 10~^. This circumstance is in common with the mo- 
mentum distribution function, kinetic energy and inter- 
nal energy. We have confirmed that 2048 Matsubara fre- 
quency is sufficient for the following results. We show the 
results obtained by using 64 x 64 meshes or 128 x 128 
meshes in the first Brillouin zone in the hole-doped case. 
We have confirmed that the finite size effects are negli- 
gible in these calculations. In the electron-doped region, 
we use 256 x 256 meshes in order to ensure the numerical 
accuracy. 

3.1 Kinetic energy 

First, Fig. 2(a) show the typical [/-dependence of SEy, 
in the hole-doped region. We have also performed the 
weak coupling theory using the second order pertur- 
bation theory (SOP) and random phase approximation 
(RPA). While the SOP and RPA are performed at T = 0, 
the FLEX is performed at T = 0.005. In the FLEX 
approximation at finite temperature, there is a criti- 
cal value of U above which superconductivity occurs. 
In the present case, Ucr/t = 2.54 at T = 0.005. This 
temperature is far below if U/t > 2.8. For example, 
Tc = 0.0102 at U/t = 4.2. 

It is clearly shown that the sign of dEk changes from 
negative to positive with increasing U/t. The negative 
sign is expected in the conventional BCS theory. In the 
BCS theory, the kinetic energy increases owing to the 
particle-hole mixing which is essential for the Cooper 
pairing. This increase of kinetic energy is logarithmically 
divergent for the cut-off of energy as 5E]^ = —pb? log 
whose absolute value is miich larger than the condensa- 
tion energy = \pl^ . Here, p is the electronic DOS at 
the Fermi level. In the SOP and RPA, these weak cou- 
pling behaviors are reproduced. We see that the FLEX 
approximation also reproduces the weak coupling behav- 
iors in the weak coupling region. 

On the other hand, the positive sign in the strong cou- 



(a) 




(b) 

— 1 — 




Fig. 2. The difTerciicc of kinetic energy between the normal state 
and the superconducting state, (a) The results of SOP, RPA and 
FLEX at <5 = 0.1. (b) The results for 10%, 15% and 20% hole- 
dopings as well as 10% electron-doping. 



pling region is a remarkably unconventional. In order to 
clarify the microscopic origin of this behavior, we show 
the momentum distribution function in Fig. 3. We see 
that n^(fc) — n^ik) takes large absolute value around the 
Fermi surface. This is because the quasi-particles near 
the Fermi surface mainly contribute to the supercon- 
ductivity. The qualitatively different behavior of Sn{k) 
in the direction perpendicular to the Fermi surface is a 
key to understand the unusual behavior. It is shown that 
n^(fc) — n^{k) is positive (negative) below (above) Fermi 
surface at the cold spot around k = (7r/2, 7r/2). The situ- 
ation is opposite at the hot spot around k = (tt, 0). This 
result means that the kinetic energy arising from the hot 
spot increases owing to the particle-hole mixing induced 
by the superconducting gap. On the other hand, the ki- 
netic energy arising from the cold spot decreases owing 
to the feedback effect on the spin fluctuation.'''* The gap 
in the magnetic excitation induced by the superconduct- 
ing gap remarkably decreases the correlation effects on 
the low energy electron states. Therefore, quasi-particles 
recover their coherent character below Tc. Since the su- 
perconducting gap is small at the cold spot, this feed- 
back effect dominates the role of particle-hole mixing. 
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The sign of SEy^ is determined by these two competing ef- 
fects. In the weak coupUng region, the particle-hole mix- 
ing is dominant and 5E\^ is negative. On the other hand, 
the feedback effect is dominant in the strong coupUng 
region where 5Ei^ is positive. The qualitatively similar 
effect has been discussed phenomenologically as a "quasi- 
particle undressing". ^^'^^'^^'^^'^^ Here, the "undressing" 
is caused by the feedback effect on the spin fluctuation. 
Note that the positive sign of 5E]^ has been reported in 
the spin-fermion modc^r^^^*" where the momentum de- 
pendence is simply neglected. However, our microscopic 
calculation shows that the momentum dependence plays 
an essential role. Note that the positive value of 6E]^ is 
also consistent with some numerical methods. ^^'^^ 

The positive sign of SEy^ means that the kinetic en- 
ergy plays a role for lowering the internal energy be- 
low Tc- This is in sharp contrast to the BCS the- 
ory, but as shown here, the spin-fluctuation-induced- 
superconductivity also gives kinetic energy gain. It 
should be noticed that the FLEX approximation is a 
microscopic description of the nearly anti-ferromagnetic 
Fermi liquid theory, where the Cooper pairing between 
quasi-particles occurs. Therefore, the kinetic energy gain 
is not a consequence of the non- Fermi liquid normal state 
as assumed in Ref. 39, and it is not a negative evidence 
for the concept of Cooper pairing between quasi-particles 
as argued in Refs. 28-30. 

Generally speaking, the kinetic energy gain can occur 
in the strong coupling superconductors where the feed- 
back effect is important. However, such a strong feedback 
effect as to change the sign of dEy is not expected in the 
low-Tc superconductors. We point out that the gain of 
kinetic energy in Fig. 2 is due to (i) strong AF spin fluc- 
tuation, (ii) high-Tc, namely large superconducting gap 
and (iii) d-wave symmetry. The existence of line node 
due to (iii) is especially important as is shown in Fig. 3. 

Next, we discuss the doping dependence. Fig. 2(b) 
shows the results of SE^ for various dopings. We have 
shown the results in the electron-doped region as well as 
in the hole-doped region. Here, the temperature is fixed 
to be T = 0.005 in the hole-doped case and T = 0.003 in 
the electron-doped case, respectively. We see that the ki- 
netic energy gain occm's in the hole-doped region and the 
crossover value of U increases with increasing the doping. 
Thus, the kinetic energy gain is likely in the under-doped 
region. 

It is interesting that the kinetic energy gain does not 
occur in the electron-doped region. The region of super- 
conducting state is very narrow as U/t = 3.8 ^ 4, and the 
tendency to the superconductivity is remarkably weak in 
the electron-doped region. This is mainly because the 
electronic DOS is small and also because the spin fluctu- 
ation is sharply localized in the momentum space. ^^'^^ 
The latter indicates that the spin fluctuation is very 
weak. Actually, an anti-ferromagnetic instability occurs 
&t U/t > 4 if we introduce a weak three-dimensionality. 
The absence of kinetic energy gain is due to the sig- 
nificantly small magnitude of superconducting gap. If 
we define A as the maximum of superconducting gap 
A(fe), we obtain nearly the BCS value 2A/Tc = 4 ~ 5 
in the electron-doped region, while 2A/Tc = 8 ~ 10 in 
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Fig. 3. Contour plot of the difference of momentum distribution 
function between the superconducting state and normal state, 
n^(fc) - nN(fc). We show the result for U/t = 4.2, S = 0.1 and 
T = 0.005. 



the under-doped region. Combined with small value of 
Tc, the superconducting gap A in the electron-doped re- 
gion is much smaller than that in the hole-doped region. 
Therefore, the effect of superconducting gap on the spin 
fluctuation is not so significant. We find that neither the 
magnetic resonance peak nor the kink in the electronic 
dispersion, which are interesting subsequences of feed- 
back effect, appear in the electron-doped region. 

3.2 Optical sum rule 

The role of kinetic energy has been discussed exten- 
sively with the relation to the optical sum rule which is 
measured experimentally.^-'^' ^^'^""^^ Here, we show that 
the decrease of kinetic energy can be observed by the 
measurement of optical integral, although they are dif- 
ferent quantities. 

In the isotropic system Hke '^He, the optical integral, 
namely the frequency integral of optical spectrum, is con- 
served through the superconducting transition as. 



f 

J — 1 



CTxx(w)da; = we — . 

m 



(26) 



This is called Ferrell-Grover-Timkam sum rule.^^ Here, 
(Txx(w) is the optical conductivity including the 6- 
function at w = 0, which corresponds to the superfluid 
density. This sum rule is generally violated under the 
periodic potential, namely in the metals. According to 
the Kubo formula,®'* the optical integral is related to the 
momentum distribution function n{k) as,^^ 



crxx(t^)da; = 27re Efe- 



-n(fc). 



(27) 



If we assume t' = 0, eq. (27) is expressed by the kinetic 
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energy as, 



r 

o — c 



(Txx(w)dw = TT^V.-, 



(28) 



This relation enables the kinetic energy to be measured 
experimentally. 

However, the optical integral is not expressed by the 
kinetic energy in more general case t' ^ 0. Therefore, 

we define the optical energy as i?op = — 2Sfe(^^^^ + 
^ dk^^ )nk so that the optical integral is expressed as, 

=2 



/ 



we 



crxx(w)da; = — T^Kp- 



(29) 



Note that E^p is equivalent to E]^ when t' = 0. 




Fig. 4. (7-dependcncc of the optical energy 5Eop for 5 = 0.1 and 
T = 0.005. We show iSEk for a comparison. 



We show the f -dependence of SEop in Fig. 4. It is 
clearly shown that the qualitative behavior of SEop is 
the same as that of (5£^k for a reaUstic value of t' /t. The 
crossover of the sign occurs in common. We see that the 
absolute value of 6Eop is larger than SE]^. However, this 
tendency depends on the value of long range hoppings.^^ 

Recently, a violation of FGT sum rule has been 
actually observed in measurements of optical inte- 
gral. ^""^'^^'^"'^-^ The results seem to be controversial, but 
some of them indicate a decrease of kinetic energy.^^'^^ It 
seems that a very accurate measurement is needed since 
the change of optical integral may be much smaller than 
its absolute value. 

Note that the maximum value of SEop shown in Fig. 4 
is 6Eop ~ 0.2meV if we adopt the band width W = 
8t = 2eV. This value is smaller than the reported value, 
6Eop ~ ImeV.^^ However, this can be a sufficient agree- 
ment because the experimental value significantly de- 
pends on the frequency cut-off^^ which is needed for a 
validity of the single-band description. 

3.3 Free energy and internal energy 

At the last of this section, we analyze the 
thermodynamic properties of spin-fluctuation-induced- 
superconductivity in more details. First, we show the 



temperature dependence of internal energy SE, kinetic 
energy SEk, optical integral SEop as well as the free 
energy SF in Fig. 5. We estimate the free energy by 
eq. (9) and F = + jin. Since the FFT involves un- 
physical results around the cut-off frequency, we replace 
the summation of Matsubara frequency in eq. (9) as 
X!-JVf% ^ YZ'n^h- This ingenuity significantly im- 
proves an accuracy of numerical calculation. Here, we 
estimate the entropy S = hy polynomial fitting and 
obtain the internal energy as £^ = F -|- TS. We show the 
result for 5E only around T = Tc because the polynomial 
fitting is not so accurate at low temperatures. 

As is shown in the figure, we obtain the condensation 
energy SF ^ 5 x 10"'' ~ 0.25meV which is consistent 
with experimental value. As for temperature depen- 
dence, the gain of kinetic energy shows qualitatively sim- 
ilar behavior to that of internal energy. The decrease of 
SE]^ as decreasing temperature is mainly owing to the de- 
crease of kinetic energy in the normal state. The damp- 
ing of quasi-particle in the normal state is reduced by 
decreasing temperature, especially at the cold spot. This 
temperature dependence of SE]^ is qualitatively indepen- 
dent of the parameters. Although SE^^ at T = seems 
to be negative in Fig. 5, we find that SE]^ at T = is 
positive when U /t and t'/t are large. 

It is shown that a considerable part of the gain of in- 
ternal energy is attributed to the kinetic energy, espe- 
cially around T = T^. However, we see (5£^k < SE in 
the whole temperature region. This means that the cor- 
relation energy also plays a positive role for the gain of 
internal energy. These features are robust in the interme- 
diate coupling region U/t = 3^5, but the contribution 
of kinetic energy increases with increasing U jt. We find 
that most part of the condensation energy is attributed 
to the kinetic energy at i' jt = 0.35, U/t = 7 and T = 0, 
but the kinetic energy gain is still smaller than the con- 
densation energy. In the variational Monte Carlo study 
for the Hubbard model, there is a parameter region 
10 < U/t < 12 where the kinetic energy and correlation 
energy play cooperative role. 




LJJ 

to 2x10 



0.002 0.004 0.006 0.008 O.C 

T 



Fig. 5. Temperature dependence of SE, 5Ek, ^^op and SF for 
5 = 0.1 and CZ/t = 4.2. 
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It should be noticed that the temperature dependence 
of optical integral 6Eop is different from SEi^, qualita- 
tively. Thus, it is necessary to distinguish the optical 
integral and kinetic energy when detailed properties of 
optical integral are discussed experimentally. As we have 
shown in Fig. 3, SE]^ and dEop are determined by com- 
petitive contributions from the hot spot and cold spot. 
Therefore, the detailed properties are sensitive to the 
long range hoppings. 

The specific heat over tempcratm-c is estimated from 
the free energy as C/T = — and wc obtain a very 
large jump of specific heat as {C^ —C^)/C^ 5.3, while 
(C^ - C^)/C^ ~ 1 in the weak coupling rf-wave BCS the- 
ory. This enhancement of specific heat jump is basically 
caused by the rapid increase of superconducting gap be- 
low Tc, which is a characteristic property of strong cou- 
pling superconductors. It should be noted that the FLEX 
approximation is not appropriate in the under-doped re- 
gion around T =Tc, because the superconducting fluc- 
tuation plays an important rolc.^^ Therefore, the large 
jump of specific heat is not observed in the under-doped 
region. However, this large jump has been observed in 
optimally-doped region^'' as well as in a heavy fermion 
resemblance CeMIns.^* where the pseudogap phenomena 
hardly occur. 

Here, we propose another interpretation of the con- 
densation energy. According to eq. (9), the condensation 
energy is expressed as 6F = SV,b + SQp + SQo + nSfi. 
Fig. 6 shows the contributions from the first and second 
terms. The first term is expressed by eqs. (11) and (19). 
If we ignore the contribution from the charge susceptibil- 
ity, eq. (19) is equivalent to the free energy arising from 
the spin fluctuation discussed by Brinkman et al.^^ They 
have estimated this term within RPA and shown that the 
Anderson-Brinkman-Morel (ABM) state in ^He is stabi- 
lized by the feedback effect. The expression of eq. (19) is 
also equivalent to the free energy discussed in the SCR.^° 
Therefore, the first term S^b can be regarded as a free 
energy arising from the spin fluctuation. Fig. 6 show that 
the first term 5^1b positively contributes to the conden- 
sation energy, while the second term Silp is negative at 
T > 0.003. We see that the magnitude of 6Qf is very 
small at low temperatures and therefore the condensa- 
tion energy is basically determined by the contribution 
from SflB-"^^ This result implies that the condensation 
energy mainly originates from the feedback effect on the 
spin fluctuation. Note that the large spin fluctuation gen- 
erally lowers Qb . Although the static spin susceptibility 
Xsiq) at ifln = is reduced by the superconducting gap, 
the dynamical part at if2„ is enhanced by the feed- 
back effect. In the present case, the contribution from 
the dynamical part over-compensates the static part and 
induces the condensation energy of superconductivity. 

Note that this interpretation is simple, but not unique, 
since the classification of free energy in eq. (9) is some- 
what arbitrary. The spin fluctuation also affects the sec- 
ond term Sflp through the self-energy, and this contribu- 
tion is necessary so as to satisfy the conservation laws. 
Nevertheless, this interpretation of condensation energy 
is somewhat interesting. We see in eq. (22) that the ef- 



T ' 1 ' r 




Fig. 6. Temperature dependence of condensation energy 5F, and 
contributions from SQp and SQb- We choose 5 = 0.1 and U/t = 
4.2. 



fcctivc interaction mediated by the spin fluctuation in- 
duces the d-wave superconductivity. Then, the static part 
of spin fluctuation mainly works as a de-pairing effect 
through the normal self-energy, while the dynamical part 
works as a pairing effect through the effective interac- 
tion. This frequency dependence is qualitatively consis- 
tent with the above interpretation on the condensation 
energy. 

We think that this interpretation of condensation en- 
ergy is qualitatively similar to the previous proposal^^^^* 
on the anti-ferromagnetic exchange energy arising from 
the magnetic resonance peak. In the FLEX approxi- 
mation, the magnetic resonance peak clearly appears. -'^^ 
Then, the frequency dependence of spin susceptibility on 
the real axis is determined by the frequency dependence 
on the imaginary axis. The latter has been discussed 
above and then we found that the dynamical part with 
if2„ ^ increases owing to the superconductivity. This 
increase is an origin of magnetic resonance peak appear- 
ing in the spin susceptibility on the real axis. There- 
fore, the interpretation of condensation energy discussed 
above is directly related to the appearance of the mag- 
netic resonance peak. 

At the last of this section, we discuss the doping de- 
pendence of condensation energy. As shown in Fig. 7, the 
condensation energy has a dome shape. This behavior 
should be contrasted to the fact that increases with 
under-doping. The decrease of condensation energy in 
the under-doped region is basically due to the existence 
of competing order, namely the anti-ferromagnetism in 
the present case. In the under-doped region, the static 
spin correlation remarkably decreases the free energy in 
the normal state F^, and therefore the condensation en- 
ergy SF = F^ — F^ is reduced. As shown in the inset. 
6F increases monotonically as a function of U. Although 
the phenomenological treatment has concluded that the 
dome shape appears by increasing the effective coupling 
constant,^" our microscopic theory provides a qualita- 
tively different result. We have confirmed that these pa- 
rameter dependences are consistent with the numerical 
estimation of eqs. (17) and (18). For example, /x^ — /U^ 
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in eq. (17) changes its sign around S = 0.125. 



3x10 



O 
5J- 



2x10 



o 



CO 1x10 




0.25 



Fig. 7. Doping dependence of condensation energy SF for U/t = 
4.2 and T = 0.005. Tc/iO is shown for a compajrison. The inset 
shows the {/-dependence of condensation energy for S = 0.1 and 
T = 0.005. 



The dome shape of the condensation energy is quali- 
tatively consistent with experimental observations.®^' 
However, we wish to stress that our definition of conden- 
sation energy is somewhat different from that used in the 
experimental analysis. Since a pseudogap exists above Tc, 
the extrapolation of the normal state free energy to T < 
Tc needs many cares. If the superconducting fluctuation 
is an origin of pseudogap phenomena, '"''^ the free energy 
is reduced by the fluctuation above T^. In this case, the 
condensation energy may be under-estimated if the role 
of pseudogap is neglected.''^ Actually, the experiments 
have observed a remarkable suppression of condensation 
energy in the under-doped region.®®'®^ 

4. Summary and Discussions 

In this paper, we have analyzed a mechanism of high- 
Tc siiperconductivity with a main interest on the ener- 
getics. The Hubbard Hamiltonian was analyzed on the 
basis of the FLEX approximation. It is shown that the 
kinetic energy is decreased below Tc in the under-doped 
region, while it is increased in the over-doped region. In- 
terestingly, the gain of kinetic energy can not occur in 
the electron-doped region. 

These findings are related to the violation of the opti- 
cal sum rule. The recent observation in the under-doped 
region has reported incompatible results to the BCS the- 
ory. However, as shown here, the microscopic theory 
on the spin-fluctuation-induced-superconductivity repro- 
duces the decrease of kinetic energy in relatively strong 
coupling region, U/t > 4. The decrease of kinetic energy 
in the ordered state is quite unusual, because the phase 
transition from the normal state is usually induced by 
the correlation energy which competes with the kinetic 
energy. The microscopic origin of the decrease of kinetic 
energy is the feedback effect.^' ^ The low energy spin fluc- 
tuation significantly suppresses the coherent motion of 
quasi-particles above Tc, and therefore quasi-particle life- 
time is short. Below Tc, the low energy spin fluctuation 



is suppressed by the opening of superconducting gap. 
Therefore, the coherence of quasi-particles is recovered 
below Tc especially around (7r/2,7r/2), and the kinetic 
energy decreases. 

We have shown that another interpretation of conden- 
sation energy is possible. As discussed in §3.3, the con- 
densation energy is dominated by SflB which is expressed 
as eq. (19). This term has been discussed as a contribu- 
tion from the spin fluctuation to the condensation en- 
ergy.®^ We flnd that while the static part of spin fluctu- 
ation plays a negative role for the condensation energy, 
the dynamical part plays a positive role. This interpre- 
tation is complementary with the analysis of the kinetic 
energy discussed above. 

Although we have discussed the mechanism of high- 
Tc superconductivity from the viewpoint of the conden- 
sation energy, we wish to stress that the understanding 
obtained from the analysis of effective interaction is clear 
and useful. As shown in this paper, the microscopic origin 
of the energy gain depends on the doping or U/t. For ex- 
ample, (5i?k is positive in the under-doped region, while it 
is negative in the over-doped region. On the other hand, 
the dominant scattering process leading to the supercon- 
ductivity is universal, namely the anti-ferromagnetic spin 
fluctuation induces the strong scattering from (±7r, 0) to 
(0, ±7r) which is attractive in the d-wave channel. In other 
words, even if the behaviors of kinetic energy are differ- 
ent between under-doped and over-doped region, it does 
not mean that pairing mechanism changes as a function 
of doping. 
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Appendix: Derivation of eqs. (9-11) 

Here, we provide a derivation of eqs. (9-11) which are 
the general relation between the thermodynamic poten- 
tial and Green function in the superconducting state. 
The derivation in the normal state has been given by 
Luttinger and Ward."'^ This can be extended to the su- 
perconducting state in a straightforward way. 

We adopt the "Bogolyubov's trick" in order to for- 
mulate the perturbation theory in the symmetry broken 
state. We add an infinitesimal symmetry-breaking term 
in the Hamiltonian as H —y H' , 



Hi, 



(A-l) 



K = J2<''K.''ka - ^o^T^^-fei - ^Sc_fe^Cfe^,(A.2) 



Hi = U ^Ui-iriii. 



(A-3) 



Here, Aq is an infinitesimal value which breaks the U{1) 
gauge symmetry. We take the limit Aq — > at the last 
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of the derivation. 

We denote the Green functions and self-energy in a 
matrix form. 



G{k) 



S(fc) 



G{k) F{k) 
F^k) -G{-k) J ' 



(A-4) 



(A-5) 



-A*(fc) -E"(-fc) 
Here, the non-interacting Green function is obtained as 

G(°)(A:)-i Ao 



Go(fc) 



(A-6) 



A* _G(o)(-fc)-i 
Then, Dyson-Gorkov equation is described as, 

G{k)^Go{k) + Goik)±{k)G{k). (A-7) 

Following the perturbation theory with respect to Hj, we 
obtain the differential of thermodynamic potential with 
respect to t/, 



(A-8) 



where Si (fc) is the reducible self-energy. According to the 
equation T,i{k) = E(fc) -|- T,{k)Go{k)'Ei{k), we obtain, 

k 

On the other hand, the differential of the generating func- 
tion $[G] defined by the skeleton diagram is obtained as. 



_ J_ 
dU ~ 2U 



TrG{k)t{k) + J2 ^ . (A-11) 



2U 

k k 

Following oqs. (A-9) and (A-11), 

k 



(A-12) 



.$:TV(Go(fc)-^-G(fc)-^)^. (A-13) 



Integrating the right hand side of eq. (A-13) with respect 
to U, we obtain 

O _ $ _ n„ = _ ^[iog(i£^g(^) + TrG(A:)E(A:)]. 



detGo(A;)-i' 



(A-14) 



Here, we have used the relations f2 = fio and $ = at 
f/ = 0. Taking the limit Aq — > 0, we obtain eqs. (9-11). 
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